clc,clear
f = @(x,y) x.^2 + 5*sin(y);
g = @(x,y)[2*x;
        5*cos(y)];
[x, val, k] = FR_CG(f, g, [5; 5]);
x
val
k